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Abstract 


We  develop  a numerical  method  to  study  the  motion  of  a planar  crystal-melt  in- 
terface during  the  directional  solidification  of  a binary  alloy  in  the  presence  of  a time- 
dependent  electric  current.  The  model  includes  the  Thomson  effect,  the  Peltier  efiFect, 
Joule  heating  and  electromigration  of  solute  in  the  coupled  set  of  equations  governing 
heat  flow  in  the  crystal  and  melt,  and  solute  diffusion  in  the  melt.  For  a variety  of 
time- dependent  currents,  the  temperature  fields  and  the  interface  velocity  are  calcu- 
lated as  functions  of  time  for  indium  antimonide  and  bismuth,  and  for  the  binary  alloys, 
germanium-gallium  and  tin-bismuth.  For  the  alloys,  we  also  calculate  the  solid  com- 
position as  a function  of  position,  and  thus  make  quantitative  predictions  of  the  effect 
of  an  electrical  pulse  on  the  solute  distribution  in  the  sohdified  material.  In  addition, 
for  a sinusoidal  current  of  small  amplitude,  we  have  compared  the  numerical  solutions 
with  approximate  analytical  solutions  valid  to  first  order  in  the  current  amplitude.  By 
using  this  numerical  approach  the  specific  mechanisms  which  play  dominant  roles  in 
interface  demarcation  by  current  pulsing  can  be  identified. 
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1.  Introduction 


A number  of  effects  occur  when  an  electric  current  is  passed  through  a solidifying  sample.  In 
the  bulk  of  the  sample,  there  is  Joule  heating  and,  in  the  presence  of  a temperature  gradient, 
there  is  Thomson  heating  or  cooling.  In  addition,  at  the  crystal- melt  interface  there  is  lib- 
eration or  absorption  of  heat  due  to  the  Peltier  effect.  The  Thomson  and  Peltier  effects  are 
proportional  to  the  current  and  thus  depend  on  its  sign,  while  Joule  heating  is  proportional 
to  the  square  of  the^jb^rrent.  In  alloys,  electric  currents  also  cause  electromigration  which 
modifies  the  solute  distribution.  Thus,  the  passage  of  an  electric  current  through  the  sample 
during  crystal  growth  alters  the  compositional  and  thermal  profiles  in  the  crystal  and  in 
the  melt,  and  also  modifies  the  interface  velocity,  thereby  changing  the  final  properties  of 
the  solidified  sample.  We  numerically  model  the  effect  of  a time-dependent  electric  current 
on  the  directional  sofidification  of  a binary  alloy  with  a planar  crystal-melt  interface.  The 
calculations  elucidate  the  relative  importance  of  various  electrical  phenomena  on  the  solid- 
ification process,  namely,  the  roles  of  the  Thomson  effect,  the  Peltier  effect.  Joule  heating 
and  electromigration  of  solute. 

There  have  been  a number  of  studies  of  the  effect  of  an  electric  current  on  solidification 
and  crystal  growth.  Early  work  on  the  effects  of  an  electric  current  on  solidification  has  been 
described  by  Pfann[l].  The  Peltier  effect  between  crystal  and  melt  was  measured  in  bismuth, 
cadmium  and  zinc  by  passage  of  a direct  current  by  Bardeen  and  Chandrasekhar [2].  The 
effect  of  a direct  current  on  solute  redistribution  in  bismuth- tin  alloys  and  on  interfacial 
instability  in  tin-bismuth  alloys  has  been  studied  by  Verhoeven[3,  4].  More  recently,  electric 
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currents  have  been  used  for  the  growth  of  epitaxial  layers  of  gallium- arsenide  (liquid  phase 
electroepitaxy,  see,  for  example,  reference [5]).  The  elTect  of  an  electrical  current  on  morpho- 
logical stability  has  been  the  subject  of  a number  of  recent  theoretical  investigations [6,  7,  8,  9]. 

Singh  et  ah  [10]  employed  current  pulsing  to  mark  the  crystal-melt  interface  in  tellurium- 
doped  indium  antimonide  and  thus  to  measure  crystal  growth  rates.  This  technique  has 
been  used  to  study  interface  instability  in  gallium-doped  germanium[ll]  and  to  measure 
the  Peltier  coefficient  for  indium  antimonide[12].  Peltier  heating  due  to  a current  pulse 
in  a pure  material  has  been  numerically  modelled  by  Wargo  and  Witt  [13];  this  model  has 
recently  been  extended  by  Duffar  et  ah  [14]  Thermal  measurements  in  pure  bismuth  and 
bismuth-manganese  eutectics  during  electrical  pulsing  have  been  carried  out  by  Silberstein 
et  aL[15,  16]. 

The  equations  governing  alloy  solidification  are  inherently  nonlinear  since  the  crystal- 
melt  interface  is  a free  boundary  (see,  for  example,  Derby  and  Brown[17]  and  references 
therein).  We  numerically  treat  the  effect  of  a time  dependent  electric  current  on  the  planar 
directional  solidification  of  a binary  alloy  by  performing  a transformation  of  variables  which 
fixes  the  solid  and  fiquid  regions  in  the  transformed  domain.  The  spatial  derivatives  of  the 
transformed  differential  equations  are  replaced  with  second  order  accurate  finite  differences, 
and  the  resulting  set  of  equations  are  then  integrated  in  time  using  ordinary  differential 
equation  integration  software.  For  a variety  of  time-dependent  currents,  the  temperature 
fields  and  the  interface  velocity  as  functions  of  time  are  presented  for  indium  antimonide 
and  bismuth,  and  for  the  binary  alloys,  germanium-gallium  and  tin-bismuth.  For  the  alloys. 
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we  also  calculate  the  solid  composition  as  a function  of  position,  and  thus  make  quantitative 
prediction  of  the  effect  of  an  electrical  pulse  on  the  solute  distribution  in  the  solidified 
material.  The  model  also  allows  us  to  investigate  the  role  of  individual  thermoelectric  effects; 
results  are  presented  for  bismuth. 

2.  The  Model 

We  consider  the  one- dimensional  case  of  the  free  boundary  problem  described  by  Coriell 
et  al.[9],  in  which  the  electric  current  density  I(t')  is  allowed  to  vary  with  time,  t',  and  the 
crystal-melt  interface  moves  along  the  z'  axis.  The  governing  set  of  differential  equations  for 
the  temperature  T(z\t')  and  for  the  composition  C(z',t')  in  the  bulk  Hquid  phase  are 

OTl  . d^Ti  , P ^OTl 

dt'  cTi  dz'  ’ 

dCi  d'^Ci  I dCi 

dt'  ~ ^ dz'^  “cTi,  dz'  ’ 

while  in  the  solid  phase 

dTs  , d'^Ts  P ^dTs 

governs  heat  flow.  Diffusion  of  solute  in  the  sohd  phase,  the  density  change  upon  soHdiflcation 
and  fluid  flow  will  be  neglected.  In  Eqns.  (1)  to  (3),  Cp  is  the  speciflc  heat  per  unit  volume, 
A is  the  thermal  conductivity,  cr  is  the  electrical  conductivity,  r is  the  Thomson  coefficient 
and  u is  the  differential  mobihty  of  solute,  with  subscripts  5 and  L denoting  sohd  and  hquid, 
respectively. 

The  compositional  and  the  thermal  fields  are  coupled  at  the  crystal- melt  interface  z'  = 
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h(t')^  through  the  boundary  conditions 


Ts-Tl  = 0, 


(4) 


Tl  = Tm  + rnCi^ 


m 

dz' 


Dl 


, fdTA  dh  ^ ^ 


dCi 


uCl  ^ dh 


(5) 

(6) 
(7) 


^ dz'  J (7l  dt' 

in  which  k = CsjCi  is  the  solute  distribution  coefficient,  Tm  is  the  melting  point  of  the  pure 
material,  m is  the  slope  of  the  hquidus,  IIsl  is  the  Peltier  coefficient,  and  Zt,  is  the  latent 
heat  of  fusion  per  unit  volume.  All  of  the  thermophysical  properties  in  eq.  (1)  through  eq. 
(7)  are  assumed  constant. 

During  the  sohdification  of  a binary  alloy,  the  interfacial  temperature  is  a priori  unknown, 
and  must  be  determined  self-consistently  from  the  phase  diagram  as  given  by  eq.  (5).  We 
note  that  to  accurately  model  the  effects  of  electrical  currents  sufficiently  large  to  cause 
melting  of  an  alloy,  diffusion  in  the  crystal  must  be  considered.  Numerical  solution  to  the 
set  of  equations  (1)  to  (7),  together  with  appropriate  far-field  boundary  conditions  (discussed 
below),  will  determine  the  position  of  the  crystal- melt  interface  as  a function  of  time,  and 
will  provide  the  fully  time-dependent  behavior  of  the  thermal  and  compositional  profiles 


throughout  the  sample. 
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3.  Numerical  Method 


3.1.  Variable  Transformations 


Our  numerical  solution  procedure  is  implemented  by  transforming  the  set  of  equations 
(1)  to  (7)  to  a moving  coordinate  system  traveling  with  the  sample  at  velocity  V,  such  that 
2:  = z'  — Vt'  and  t = t' . In  the  moving  coordinate  system  the  interface  position  w[t)  is 
given  by  w[t)  = h{t')  — Vt' . The  extent  of  both  the  liquid  and  solid  regions  is  truncated 
by  introducing  finite  cut-off  lengths,  where  far-field  boundary  conditions  are  appHed.  We 
have  chosen  the  origin  of  our  coordinate  system  such  that  the  solid  phase  lies  within  the 
region  0 < 2:  < w{t).  We  define  a thermal  cut-off  length  2 = Lt  and  a solutal  cut-off 
length  ^ = Lc.  Since  the  thermal  and  solutal  fields  are  only  coupled  at  the  crystal- melt 
interface,  it  is  not  necessary  to  take  Lt  = Lc\  in  order  to  accurately  resolve  the  solute  field 
and  save  computational  time,  Lc  will  be  chosen  smaller  than  Tj.  We  fix  the  position  of  the 
interface  by  mapping  the  crystal  and  the  melt  into  unit  intervals  using  the  transformations 
= zlw[t)  in  the  solid,  = [z  — w[t)]l[LT  — n;(t)]  + 1 in  the  liquid  for  the  thermal  field 
and  ic  = [z  — n^(^)]/[^c  “■  n;(t)]  + 1 in  the  liquid  for  the  compositional  field. 

The  equations  governing  heat  flow  and  diffusion  in  the  melt  become, 
dTi  \l 


dt  Cpl{Lt  - wy  dil 
dCi 


dTL\  V-u;,(a-2)-ri//cpi  /arA  P 

\diLr 


(Lt  — w) 


O'LCpL 


^ Dl  (d^CA  V - w,(A  - 2)  - uIAt  (dCA 

dt  {Lc-wy[dec)  (Lc-w)  [die)’ 

respectively,  and  in  the  solid  the  equation  governing  heat  flow  becomes 

dTs 


(8) 

(9) 


dt 


As  fd^Ts\  ^ V + w,is  - rsl/c^s  fdTs\  P 


CpSW^  V 


w 


O' SCpS 


(10) 


with  Wt  = dw/dt.  At  the  crystal-melt  interface  (<^5  = = 1),  eqs.  (4)  and  (5)  remain 

unchanged,  and  eqs.  (6)  and  (7)  become: 


A. 


dT, 


A. 


(dTi 


LyW  \ d^s  J 

Dl 


L.[Lt 



(Lc  ~ w)  \ d^c  , 


uCl 


, , Hsl/ 

Wt  -\-  V H — , and 


(^L 


I-{V  + Wt)(l-k)CL. 


(11) 

(12) 


Finally,  the  followdng  mixed  boundary  conditions  are  applied  at  the  far- field  cutoff  lengths 


400  / ^'T  \ 

(^y)+Brr.  = Erata  = 2,  (13) 

Aoo  / QQ  \ 

(17^  j ^ 

= = (15) 

where  the  quantities  with  superscript  ’00’  are  constants. 


3.2.  Discretization 

We  reduce  the  coupled  set  of  partial  differential  equations  and  interfacial  conditions  to  a 
system  of  ordinary  differential  equations  in  time  by  discretizing  the  spatial  derivatives.  The 
spatial  derivatives  in  each  of  the  differential  equations  (8)  to  (10)  are  replaced  by  second 
order  accurate  centered  finite  differences  at  equally  spaced  interior  nodes.  At  the  exterior 
boundaries  we  apply  one-sided  second  order  accurate  finite  differences  as  approximations 
to  the  first  order  derivatives  in  eqs.  (13)  to  (15),  and  eliminate  the  temperature  and  the 
composition  as  unknowm  variables  at  the  sample  endpoints. 

The  conservation  of  solute  boundary  condition,  eq.  (12),  is  treated  by  introducing  a 
fictitious  nodal  point  in  the  solid  phase  and  then  using  a centered  difference  approximation 
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for  the  derivative  of  melt  composition  at  the  interface.  The  fictitious  nodal  point  introduces 
an  extra  unknown  into  the  problem,  and  the  interfacial  composition  can  be  determined  by 
applying  the  diffusion  equation  for  the  bulk  liquid  phase  at  the  interface  as  well.  Finally, 
the  heat  flux  condition  eq.  (11)  provides  an  ordinary  differential  equation  for  the  interface 
position;  second  order  accurate  one-sided  differences  are  used  to  evaluate  the  temperature 
gradients  in  the  hquid  and  soHd  at  the  interface. 

Thus,  specifying  appropriate  experimental  and  material  parameters,  including  the  coeffi- 
cients in  the  far-field  conditions  (Eqns.  (13),  (14)  and  (15)),  initial  data  for  the  temperature 
and  solutal  fields  and  the  interface  position,  we  then  integrate  the  ordinary  differential  equa- 
tions in  time  using  the  software  package  SDRIV3[18].  This  determines  the  thermal  and  the 
solutal  fields  and  the  position  of  the  crystal-melt  interface  cls  functions  of  time. 

3.3.  Absorbing  Far-Field  Boundary  Conditions 

In  order  to  accurately  resolve  the  solutal  field  in  front  of  an  advancing  interface,  we  impose 
a cutoff  boundary  that  is  often  far  shorter  than  a typical  experimental  sample  length,  and  is 
only  a few  solute  decay  lengths  from  the  interface.  We  use  the  following  far-field  boundary 
condition: 

{V  - [ullai])  {w )+^^  = 

where  Coo  is  the  concentration  of  solute  in  the  liquid  far  from  the  interface.  Eq.  (16)  is 
automatically  satisfied  under  steady-state  growth  conditions  in  the  sample  at  any  value  of 
z.  Similar  boundary  conditions  could  also  be  imposed  on  the  thermal  fields  at  the  sample 
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endpoints;  however  for  the  calculations  reported  here  we  simply  fix  the  temperatures  at  the 
outer  boundaries  using  values  consistent  with  steady-state  motion. 

3.4.  Numerical  Comparisons  to  Known  Solutions 

We  compare  our  numerical  results  with  an  analytical  solution  given  by  Smith  et  al.[19] 
and  by  Tiller  and  Sekerka[20]  for  the  initial  transient  in  an  initially  homogeneous  sample 
solidifying  at  constant  velocity  with  constant  interface  temperature  (that  is,  m = 0).  In 
our  numerical  calculations  the  imposed  thermal  fields  keep  the  interface  moving  at  constant 
velocity,  and  we  calculate  the  solutal  profile  as  it  changes  from  its  initially  homogeneous  state. 
As  in  aU  of  our  calculations,  we  use  the  mixed  condition  given  by  eq.  (16)  at  the  sample 
endpoint  Lc-  We  compare  our  numerical  calculations  of  the  interfacial  melt  composition  to 
eq.  (26)  of  Smith  et  al.[19]  at  various  times  in  Table  (1)  for  different  values  of  the  mesh 
spacing.  The  numerical  and  the  analytical  results  are  in  good  agreement,  and  indicate  that 
the  numerical  method  is  spatially  second  order  accurate. 

We  also  compare  our  numerical  results  to  eq.  (43a)  of  Smith  et  al[19]  which  represents 
the  changes  that  the  solutal  field  undergoes  when  there  is  an  instantaneous  change  in  the 
interfacial  velocity  from,  say,  V\  to  V2.  The  initial  solute  distribution  corresponds  to  the 
steady-state  distribution  for  velocity  Vi.  (The  hquidus  slope  m is  stiU  assumed  to  be  zero.) 
We  calculate  the  resulting  transient  in  the  solutal  field  that  occurs  as  the  solutal  profile 
adjusts  itself  to  the  discontinuous  change  in  the  growth  speed.  Comparison  between  the 
analytically  and  numerically  calculated  interface  compositions  as  a function  of  time  is  given 
in  Table  (2)  for  Lc  = 4T)/Vi  = 0.4  cm.  Again,  the  numerical  and  the  analytical  results  are  in 
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good  agreement,  and  indicate  that  the  numerical  method  is  spatially  second  order  accurate. 
The  accuracy  of  our  numerical  technique  for  a discontinuous  change  in  the  interface  velocity 
indicates  that  we  can  accurately  calculate  the  solute  profiles  under  conditions  of  square  wave 
current  pulsing,  for  which  discontinuous  changes  in  the  growth  velocity  also  occur. 

In  the  Appendix  we  derive  a linear  theory  for  the  case  of  a small  amplitude  sinusoidal 
current  of  the  form  J = /^  -f-  in  which  uj  is  the  frequency  of  the  oscillation.  The 

solutions  are  assumed  to  be  of  the  form, 


\ 

Ts(z,t) 

' rf  (^)  ^ 

Mz^t) 

Tf\z) 

+ 

T)u\z) 

C{z,t) 

w(t)  ^ 

V ' y 

(17) 


The  unperturbed  solutions  (denoted  by  the  superscript  ’(0)’)  satisfy  the  steady-state  equa- 
tions for  a constant  current  Iq  and  velocity  V.  For  convenience  in  this  calculation  we  take  z 
= 0 to  be  the  position  of  the  unperturbed  interface;  the  quantity  S measures  the  deviation 
of  the  interface  position  due  to  the  oscillating  current.  Eqn.  (17)  is  substituted  into  the  gov- 
erning equations  and  boundary  conditions  given  by  Eqns.  (1)  through  (7)  and  the  far-held 
conditions,  Eqns.  (13)  through  (15).  We  expand  all  the  boundary  conditions  in  a Taylor’s 
series  about  the  unperturbed  interfacial  position  retaining  only  terms  that  are  linear  in  the 
coefficients  of  The  resulting  linearized  set  of  equations  is  then  solved.  A comparison 
between  analytically  and  numerically  calculated  growth  velocities  are  given  in  Table  (3), 
showing  agreement  between  the  two  calculations. 
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4.  Numerical  Results 


We  carry  out  numerical  calculations  for  the  pure  materials  indium  antimonide,  InSb,  and 
bismuth,  Bi,  and  for  the  binary  alloys,  germanium-gallium,  Ge-Ga,  and  tin-bismuth,  Sn- 
Bi.  The  material  properties  used  in  our  calculations  are  given  in  Table  (4).  For  all  our 
calculations,  the  initial  thermal  and  solutal  profiles  are  the  steady-state  profiles  given  by 
eqs.  (21)  - (23)  in  the  Appendix. 

Numerical  calculations  for  InSb  have  been  carried  out  by  Wargo  and  Witt  [13]  and  Duffar 
et  ah  [14]  The  thermal  fields  associated  with  current  pulsing  in  pure  Bi  have  been  studied 
by  Silberstein  and  Larson[16].  The  alloy  Ge-Ga  has  been  successfully  used  by  Holmes  and 
Gatos[ll]  for  interface  demarcation  by  current  pulsing.  Finally,  the  alloy  Sn-Bi  has  been  cho- 
sen for  experiments  in  microgravity  with  the  MEPHISTO  apparatus[21],  which  uses  current 
pulsing  for  interface  demarcation. 

4.1.  Indium  Antimonide 

The  plots  in  Fig.  (1)  give  the  interfacial  growth  speed  u = d/i/dt  as  a function  of  time.  In 
the  top  plot  in  Fig.  (1),  the  initial  growth  rate  is  11.33^m/s  and  the  initial  temperature 
gradient  in  the  melt  at  the  interface  is  25.0  K/cm.  A current  pulse  of  20.3  A/cm^  and 
of  one  second  duration  is  appHed  at  3.0  and  7.0  seconds.  In  our  model,  since  there  is  an 
instantaneous  jump  in  the  current  (A/),  there  is  an  instantaneous  jump  in  the  microscopic 
growth  rate  by  an  amount  Av  = —HslAI/Lv  = 14.6//m/s,  since  the  thermal  fields  cannot 
change  instantaneously.  This  discontinuity  in  the  velocity  is  captured  by  the  numerical 
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method.  As  the  thermal  field  relaxes,  the  velocity  decreases  by  0.45/im/s,  and  then  a jump 
in  the  velocity  of  the  same  magnitude  but  of  opposite  sign  occurs  when  the  current  pulse  is 
turned  off.  In  the  bottom  plot  in  Fig.  (1)  we  show  the  results  for  repeated  current  pulses 
of  9.5  A/cm^  which  are  on  for  20  s and  off  for  40  s,  for  which  the  initial  growth  rate  is 
12.0/xm/s  and  the  initial  temperature  gradient  in  the  melt  at  the  interface  is  25.0  K/cm. 
Both  calculations  shown  in  Fig.  (1)  are  in  close  agreement  with  Wargo’s  numerical  and 
experimental  results.  Furthermore,  as  calculated  by  Wargo  for  repeated  pulses,  the  maxima 
and  minima  of  the  velocity  of  the  interface  decrease  from  their  initial  values  as  the  transients 
decay;  the  average  growth  rate  is  eventually  determined  by  the  overall  pulling  rate. 

4.2.  Bismuth 

We  next  investigate  the  effect  of  an  electric  pulse  on  the  interfacial  velocity  and  the 
thermal  fields  in  the  liquid.  Furthermore,  we  isolate  the  contribution  of  the  Peltier,  Thomson 
and  Joule  effects  by  setting  the  appropriate  thermoelectric  coefficients  to  zero.  The  interface 
is  initially  static  and  the  temperatures  at  the  finite  cut-off  lengths  are  held  fixed.  The 
temperature  gradient  in  the  liquid  at  the  interface  is  initially  100  K/cm  and  the  length  Lj 
is  4 cm  with  the  interface  located  1.5  cm  from  the  solid  end  of  the  sample.  At  t = 20  s, 
a current  pulse  of  + 80  A/cm^  is  applied  to  the  sample  for  a 10  second  duration,  as  used 
in  the  experiments  of  Silberstein  and  Larson[16].  The  interfacial  velocity  and  the  change  in 
temperature  at  three  fixed  positions  in  the  liquid,  0.05,  1.25  and  2.27  cm  from  the  initial 
interface  position,  are  plotted  as  functions  of  time  in  Fig.  (2). 

Fig.  (2)  shows  that  the  temperature  change  near  the  interface  is  very  different  from  that 
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away  from  the  interface,  which  is  in  qualitative  agreement  wuth  experiment [16].  Quantita- 
tive comparison  between  the  numerical  and  experimental  results  is  not  possible  because  of 
uncertainty  in  the  physical  parameters  and  boundary  conditions  at  the  sample  endpoints. 
Figs.  (3)  and  (4)  show  the  individual  contributions  of  the  Peltier,  Thomson  and  Joule  effects 
to  the  temperature  change  and  the  interface  velocity,  respectively.  For  example,  in  the  top 
plot  for  the  Peltier  effect  in  Figs.  (3)  and  (4),  the  Thomson  coefficients  in  the  hquid  and 
the  solid  are  zero,  and  the  electrical  conductivities  in  the  solid  and  liquid  are  infinite.  As 
would  be  expected  from  the  sign  of  the  Peltier  coefficient,  a positive  current  causes  cooling 
in  the  liquid  and  an  increase  in  the  interface  velocity.  For  bismuth,  the  Thomson  coefficient 
in  the  solid  is  positive,  while  the  Thomson  coefficient  in  the  liquid  is  almost  zero  but  slightly 
negative  (estimated  from  Fig.  (6)  of  Ubbelohde[22]).  Thus  the  Thomson  effect  in  the  solid 
causes  coofing  of  the  sample,  and  the  Thomson  effect  in  the  Hquid  causes  heating.  From  the 
magnitude  of  the  coefficients,  the  cooling  in  the  sofid  dominates,  and  the  interface  velocity 
is  increased.  Thus  the  interface  moves  closer  to  the  point  at  which  the  temperatures  are 
measured  and  the  temperatures  are  lowered  even  though  there  is  a slight  Thomson  heat- 
ing in  the  liquid.  Additional  calculations  with  different  values  for  the  Thomson  coefficients 
have  verified  this  interpretation.  The  Joule  heating,  which  is  independent  of  the  sign  of  the 
current,  causes  a temperature  rise  and  a decrease  in  the  interface  velocity.  While  this  is  a 
nonlinear  problem  and  the  individual  contributions  are  not  additive,  comparison  of  Figs.  (2) 
and  (4)  indicate  that  immediately  followdng  the  application  of  the  current,  the  Peltier  effect 
dominates,  but  as  the  temperature  field  relaxes,  the  Thomson  effect  and  the  opposing  Joule 
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effect  provide  relatively  constant  contributions  to  the  interface  velocity  as  well.  Fig.  (3) 
shows  that  away  from  the  interface  Joule  heating  is  dominant,  but  the  heating  is  suppressed 
near  the  sample  end  (upper  curve  in  each  plot),  due  to  the  constant  temperature  boundary 
condition.  In  Fig.  (2),  the  change  in  temperature  decays  to  zero  after  the  current  pulse  is 
turned  off.  This  behavior  depends  on  the  boundary  conditions  at  the  sample  endpoint;  if, 
for  example,  we  had  used  no  flux  conditions,  there  would  have  been  a net  change  in  the 
temperature  of  the  sample. 

Figs.  (5)  - (7)  are  analogous  to  Figs.  (2)  - (4)  except  that  the  current  is  of  opposite 
sign.  In  this  case,  all  effects  act  to  increase  the  temperature  in  the  liquid  as  shown  in  Fig. 
(6),  with  the  largest  temperature  change  occurring  near  the  crystal- melt  interface.  At  the 

beginning  of  the  current  pulse,  the  interface  velocity  is  negative,  decreases  in  magnitude,  but 

% 

remains  negative  due  to  the  heating  until  the  current  pulse  ends. 

In  Figs.  (2)  and  (5),  the  temperature  change  near  the  interface  is  neither  symmetric  nor 
anti-symmetric  with  respect  to  a change  in  sign  of  the  current  pulse.  If  Joule  heating  were 
dominant,  symmetric  behavior  would  be  expected,  while  if  Joule  heating  were  unimportant, 
anti- symmetric  behavior  would  be  expected. 

4.3.  Germanium- Gallium 

For  the  case  of  alloy  solidification,  a current  pulse  modifies  the  compositional  profile  in  the 
liquid  as  well  as  the  interface  velocity  and  the  thermal  fields  in  the  solid  and  liquid.  The 
composition  in  the  solid  depends  on  the  composition  in  the  Hquid  at  the  interface,  and  is 
therefore  modified  by  the  current  pulse.  This  change  in  sohd  composition  is  the  basis  of  the 
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interface  demarcation  technique[10]. 

We  present  calculations  for  Ge-Ga  alloys,  under  different  processing  conditions.  First,  we 
consider  the  case  in  which  the  current  is  turned  on  and  left  on.  There  is  a transient  stage  in 
which  the  system  changes  from  its  steady  state  in  the  absence  of  a current  to  a new  steady 
state  for  a finite  current  (see  eqs.  (21)  to  (23)  of  Appendix).  The  initial  and  final  interface 
velocities  are  identical,  but  the  initial  and  final  concentration  and  temperature  fields  differ. 

In  Figs.  (8)  and  (9)  a current  of  ±25  A/cm^  is  applied  at  f = 10s.  The  bulk  alloy 
composition  is  0.1  at%  gallium,  the  initial  temperature  gradient  in  the  liquid  at  the  interface 
is  130  K/cm  and  the  initial  interface  velocity  is  38^  m/s.  The  temperature  changes  are 
measured  at  three  fixed  positions  (z  = constant)  in  the  liquid,  0.02  (lower  curves),  0.5 
(middle  curves)  and  0.91  cm  (upper  curves)  from  the  initial  interface  position  in  a moving 
reference  frame,  i.e.,  if  the  interface  continued  to  move  at  constant  velocity,  the  points  would 
remain  at  fixed  distances  from  the  interface.  For  a positive  current  the  Peltier  effect  provides 
cooling  while  the  Thompson  and  Joule  effects  provide  heating.  For  this  alloy,  the  Peltier 
effect  is  dominant  and  the  interface  velocity  increases  upon  application  of  the  current  which 
results  in  an  increase  in  the  composition  of  the  solid.  Upon  the  application  of  a negative 
current  (Fig.  (9)),  the  interface  velocity  and  the  composition  of  the  solid  initially  decrease. 

In  Figs  (10)  and  (11)  we  show  the  effect  of  different  sample  compositions  on  the  change  in 
the  temperature  in  the  liquid,  for  the  same  conditions  as  for  Figs.  (8)  and  (9),  respectively. 
In  each  figure,  the  upper  plot  is  for  pure  germanium,  and  the  middle  and  lower  plots  are 
for  germanium  containing  0.3  at%  and  0.5  at%  gallium,  respectively.  As  the  composition  is 
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increased  the  temperature  change  in  the  liquid  near  the  interface  becomes  larger,  reflecting 
the  change  in  interface  temperature  with  solute  concentration.  The  maximum  change  in 
temperature  is  attained  more  quickly  for  larger  compositions.  After  the  new  steady  state 
is  attained,  the  temperature  changes  are  different  for  positive  and  negative  currents.  For  a 
positive  current,  the  interface  attains  a new  position  which  is  closer  to  the  points  at  which  the 
temperature  change  is  calculated.  This,  together  with  the  readjustment  of  the  temperature 
profile  due  to  the  electric  current,  gives  rise  to  the  difference  in  the  temperature  changes  for 
positive  and  negative  currents. 

In  the  next  series  of  calculations,  we  study  the  effect  of  a current  pulse  in  germanium- 
gallium  with  parameters  typical  of  those  used  for  interface  demarcation[ll],  namely,  ± 80 
A/cm^  and  0.1  s in  duration.  The  initial  interface  velocity  is  38  /im/s  and  the  initial  tem- 
perature gradient  in  the  liquid  at  the  interface  is  130  K/cm.  Figs.  (12)  and  (13)  show 
the  changes  in  temperatiire  and  the  interface  velocity  as  functions  of  time,  and  the  solid 
composition  as  a function  of  distance  for  positive  and  negative  currents,  respectively.  For  a 
positive  current  the  Thomson  and  Joule  effects  provide  heating  while  the  Peltier  effect  cools 
the  sample.  From  Fig.  (12),  it  is  clear  that  the  Peltier  effect  dominates  near  the  interface, 
while  away  from  the  interface  there  is  a slight  heating.  Comparing  Figs.  (12)  and  (13),  the 
temperature  change  near  the  interface  and  the  interface  velocity  are  almost  antisymmetric 
with  respect  to  a change  in  the  sign  of  the  current,  which  is  consistent  with  the  dominance  of 
the  Peltier  effect.  The  interface  velocity  instantaneously  increases  (decreases)  upon  applica- 
tion of  a positive  (negative)  current;  it  then  decreases  (increases)  rapidly  followed  by  a slow 
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decay  (rise)  until  the  current  pulse  is  over.  The  composition  changes  monotonically  upon  the 
application  of  the  current  and  then  abruptly  begins  to  return  to  its  initial  composition  once 
the  current  is  over.  The  solute  composition  in  the  solid  has  a maximum  change  of  about 
five  percent  of  the  bulk  solute  concentration.  The  solid  composition  changes  over  a distance 
of  about  10  /xm;  there  are  particularly  sharp  concentration  gradients  associated  with  the 
beginning  and  end  of  the  current  pulse. 

4.4.  Tin-Bismuth 

In  Figs.  (14)  and  (15)  a current  of  ±50 A/ cm^  is  applied  for  5 seconds.  The  bulk  alloy 
composition  is  0.78  at%,  with  an  initial  temperature  gradient  in  the  liquid  at  the  interface 
of  1 35 A"/cm  and  an  initial  interface  velocity  of  1.67/im/s.  The  temperature  changes  are 
measured  at  three  fixed  positions  in  the  liquid,  0.03  (lower  curves),  0.75  (middle  curves) 
and  1.36cm  (upper  curves)  from  the  initial  interface  position  in  a moving  reference  frame. 
For  a positive  current  the  Peltier  effect  provides  cooling;  the  Thomson  effect  is  essentially 
negligible  for  tin.  Figs.  (14)  and  (15)  show  changes  in  the  liquid  temperature  and  the 
velocity  as  functions  of  time,  and  the  composition  in  the  solid  as  a function  of  position. 
The  similarity  in  the  thermal  profiles  for  both  positive  and  negative  currents  indicate  that 
Joule  heating  is  playing  the  important  role  in  altering  growth  conditions.  For  the  positive 
current  pulse,  after  the  initial  increase  due  to  the  Peltier  effect,  the  interfacial  velocity  sharply 
decreases  to  less  than  its  initial  value,  which  is  consistent  with  the  important  role  of  Joule 
heating.  Additional  calculations  with  Usl  = 0 show  that  the  change  in  temperature  and  the 
solid  composition  are  very  similar  to  those  shown  in  Fig.  (14).  For  negative  current  pulses. 
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the  velocity  decreases  and  remains  relatively  constant  until  the  pulse  is  terminated.  The 
composition  in  the  solid  is  changed  by  about  2 percent  of  the  bulk  value  for  the  negative 
pulse  and  somewhat  less  for  the  positive  pulse,  over  a distance  of  about  20  //m. 

In  the  calculations  the  effect  of  electromigration  has  been  included,  however,  for  the 
processing  conditions  considered  the  effect  of  electromigration  is  negligible.  When  the  di- 
mensionless quantity  \uII<jlV\  <C  1,  we  expect  that  electromigration  is  unimportant.  For 
the  parameters  used  in  the  tin-bismuth  calculations  \uII(7lV\  = 0.021  and  is  even  smaller 
for  the  germanium- gallium  calculations.  Calculations  with  u = 0 have  verified  that  electro- 
migration is  unimportant.  For  very  large  currents  and  small  velocities  electromigration  may 
be  important. 

5.  Conclusions 

We  have  considered  directional  sohdification  of  pure  materials  and  binary  alloys  with  a pla- 
nar crystal- melt  interface  in  the  presence  of  a time-dependent  electric  current.  A variety  of 
numerical  results  have  been  presented  for  growth  velocities,  thermal  fields  and  solid  compo- 
sitions. 

Our  results  for  indium  antimonide  are  in  close  agreement  with  the  numerical  results  of 
Wargo  and  Witt  [13].  For  pure  bismuth,  we  have  performed  calculations  for  conditions  similar 
to  the  experiments  of  Silberstein  and  Larson[16].  Given  the  uncertainties  in  the  physical 
properties  and  the  appropriate  far-held  boundary  conditions,  the  numerical  calculations  show 
reasonable  agreement  with  their  experimental  results.  In  our  numerical  calculations  (but  not 
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in  the  experiments),  it  is  easy  to  isolate  the  various  thermoelectric  effects  and  this  is  shown 
in  Figs.  (3),  (4).  (6)  and  (7). 

One  of  the  main  uses  of  electrical  pulsing  is  for  interface  demarcation  in  alloys.  This 
technique  is  based  on  the  change  in  solute  concentration  in  the  sohdified  material  due  to  the 
electrical  pulse.  The  calculations  for  germanium-galhum  provide  quantitative  predictions  of 
the  solute  composition  in  the  soHd  due  to  electrical  pulsing,  which  has  been  widely  used 
for  interface  demarcation.  Under  typical  experimental  conditions  for  interface  demarcation, 
the  magnitude  of  the  change  in  solid  composition  is  about  five  percent  of  the  average  soHd 
composition.  The  calculations  for  tin-bismuth  indicate  a similar  change  in  soHd  composition; 
however,  whether  such  compositional  variations  can  be  observed  depends  critically  on  the 
sensitivity  of  the  observational  technique.  Apparently  etching  is  an  extremely  sensitive 
technique  for  revealing  compositional  modulations  in  the  germanium-galhum  system. 

We  have  neglected  diffusion  in  the  solid,  which  prevents  the  consideration  of  electrical 
pulses  in  alloys  that  would  cause  melting  of  the  solid.  Since  the  thermal  and  solutal  fields  are 
only  coupled  at  the  interface  (in  the  absence  of  Soret  and  Dufour  effects),  it  is  straightforward 
to  extend  our  numerical  method  to  treat  diffusion  in  the  solid. 
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Appendix  : Linear  Theory  for  a Time  Periodic  Current 


In  a reference  system  moving  at  velocity  V,  the  set  of  differential  equations  governing 
the  unperturbed  problem  are: 
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in  which  5 = /^/((jA),  K = (cpV  — tIo)/X,  t]  = (V/D){1  — Qo)  and  Qq  = [ulo/crV].  The 
solutions  to  the  above  equations  are 
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for  the  solute,  where 
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and  Csi  is  the  composition  of  the  solid  phase  at  the  crystal- melt  interface.  The  unperturbed 


temperatures  are 
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and  Gi  is  the  unperturbed  temperature  gradient  in  the  melt  at  2r  = 0.  The  perturbed 


problem  is  governed  by  the  equations 
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The  solutions  to  the  perturbed  set  of  equations  (24)  to  (24)  are 
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for  either  of  the  thermal  fields,  and 
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The  four  interfacial  boundary  conditions,  (4)  - (7),  are  linearized,  and  yield 
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for  the  perturbed  problem.  The  far-held  conditions  for  the  perturbed  quantities  are  of  the 
form 


^oo_^  + = 0. 

oz 


(34) 


The  four  interface  conditions  plus  the  three  far-held  conditions  for  T^\  and  form 
a set  of  seven  linear  inhomogeneous  equations  for  the  seven  unknowns,  of,  of,  of, 
and  6.  These  hnear  equations  are  solved  numerically,  and  values  of  the  perturbed  quantities 
are  calculated. 
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Table  1:  Comparisons  of  the  solute  transient  given  by  Smith  et  al.[19]  with  the  numerical 
calculations  for  m = 0.  The  data  is  for  Ge-Ga  with  V = 0.019  cm/s  so  that  D/V  = .01 
cm  and  the  cutoff  length  is  given  by  Lc  = jV  = 0.04  cm.  The  initial  solutal  field  has  a 
uniform  value  of  unity.  N is  the  number  of  intervals  used  to  discretize  the  liquid  domain. 
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0.00118 

0.00029 

.10 

1.52350 

0.00485 

0.00122 

0.00030 

.12 

1.58136 

0.00501 

0.00126 

0.00032 

.14 

1.63579 

0.00518 

0.00130 

0.00033 

.16 

1.68750 

0.00535 

0.00134 

0.00034 

.18 

1.73699 

0.00553 

0.00139 

0.00035 

.20 

1.78462 

0.00570 

0.00143 

0.00036 
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Table  2:  Comparisons  of  the  solute  transient  given  by  Smith  et  ah  [19]  with  the  numerical 
calculations  for  the  case  in  which  m = 0.  The  interface  is  initially  travelling  at  the  constant 
velocity  Vi,  and  the  interface  velocity  is  instantaneously  doubled  at  t = 0.  The  interface 
composition  increases  and  then  settles  down  into  a steady  state  profile.  Comparison  of  the 
analytical  solution  is  made  with  numerical  calculations  for  different  mesh  spacings.  The  data 
is  for  an  alloy  of  Ge-Ga  wdth  an  initial  velocity  Vi  ~ 0.0019  cm/s  so  that  D/Vi  = 0.1  cm. 
The  cutoff  length  is  given  by  Lc  = /Vi  = 0.4  cm. 


t(secs) 

Analytical 

Analvtical  - Numerical 
%/ 

N = 50 

N = 100 

N = 200 

.02 

1.03284 

0.00610 

0.00141 

0.00035 

2.0 

1.08996 

0.00351 

0.00088 

0.00022 

4.0 

1.11509 

0.00351 

0.00088 

0.00022 

6.0 

1.12962 

0.00364 

0.00091 

0.00023 

8.0 

1.13880 

0.00379 

0.00095 

0.00024 

10.0 

1.14468 

0.00394 

0.00099 

0.00025 

12.0 

1.14831 

0.00408 

0.00102 

0.00025 

14.0 

1.15034 

0.00422 

0.00106 

0.00027 

16.0 

1.15116 

0.00434 

0.00109 

0.00027 

18.0 

1.15107 

0.00446 

0.00112 

0.00028 

20.0 

1.15027 

0.00458 

0.00115 

0.00029 

30.0 

1.13994 

0.00507 

0.00129 

0.00034 
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Table  3:  Comparison  of  the  analytical  amplitude  of  the  oscillatory  component  uod  of  the 
interface  velocity  with  the  maximum  and  minimum  values  of  the  numerically-determined 
oscillatory  component  of  the  interface  velocity  for  a current  of  the  form  I = /iexp(2u;t), 
with  (jj  = 27r  X 10“^.  The  analytical  results  are  calculated  from  the  linear  theory  derived 
in  the  Appendix.  The  numerical  results  are  determined  after  significant  decay  of  the  initial 
transient.  The  temperature  and  the  composition  are  fixed  at  the  outer  boundaries.  The 
thermal  sample  length  is  0.5  cm  for  both  the  solid  and  liquid,  and  the  solutal  cutoff  length 
is  a distance  DjV  away  from  the  interface  in  the  melt,  with  V = 38//m/s.  The  material 
parameters  are  for  germanium-gallium  alloys  listed  in  Table  (4). 


h 

Anal,  (mag.) 

Num.  (max) 

Num.  (min) 

0.1 

9.802  X 10-’’ 

9.993  X 10"' 

-9.994  X 10“' 

1.0 

9.802  X 10-« 

9.997  X 10-« 

-9.999  X 10-5 

10.0 

9.802  X 10-® 

9.521  X 10-5 

-9.854  X 10-5 
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Table  4:  Thermophysical  Properties 


Property 

Sn-Bi 

Ge-Ga 

Bi 

InSb 

CpL[J ! {cmh<)] 

1.0 

2.1 

1.46 

1.53 

Cps[J ! [cm^  K)\ 

1.0 

2.1 

1.46 

1.4 

Tm[K\ 

505 

1209 

549 

803 

L[J  lcm^\ 

410 

2520 

54.0 

1350 

\l[J  ! [cmK  s)] 

0.3 

0.7 

0.124 

0.13 

\s[J ! {cmK  s)] 

0.6 

0.24 

0.065 

0.046 

(Jl  [1/(^0772)] 

2.1  X 10< 

1.6  X 10“ 

7.8  X 10^ 

9.5  X 10^ 

crs[l/{Qcm)] 

4.5  X 10"* 

1.25  X 10^ 

3.8  X 10^ 

2.9  X 10^ 

rs[V7A'] 

-5.7  X 10-’’ 

0.0 

7.62  X 10“^ 

-4.0  X 10-’ 

nlv/K] 

0.0 

-2.7  X 10-“ 

-7.13  X 10-’’ 

-4.0  X 10-5 

nsL[V] 

-6.1  X 10-“ 

-9.7  X 10-2 

-0.014 

—9.7  X 10“^ 

D[cm^  / s] 

1.8  X 10-^ 

1.9  X 10-“ 

1.9  X 10-“ 

1.9  X 10-“ 

m[at.%/ K] 

-2.9 

-3.7 

k 

0.3 

0.086 

u[cm}  ! (V^)] 

-1.5  X 10-^ 

3.6  X 10-^ 
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Figure  Captions 


Figure  1 : Plot  of  the  interfacial  growth  velocity  as  a function  of  time  for  pure  InSb.  The 
upper  plot  shows  calculations  for  the  case  of  Is  current  pulses  of  20. 3 A/ cm?  applied 
at  times  3s  and  7s.  The  bottom  plot  shows  the  case  of  repeated  current  pulses  of 
9.5A/cm^  which  are  on  for  20s  and  off  for  40s.  The  initial  temperature  gradient  in  the 
melt  at  the  interface  is  25. 0/v /cm  in  both  cases. 

Figure  2 : The  change  in  temperature  and  the  interface  velocity  as  a function  of  time 
for  pure  Bi,  with  a current  pulse  of  80A/cm^  and  duration  10s.  The  temperature 
changes  are  measured  at  three  fixed  positions  in  the  hquid,  0.05  (lower  curve),  1.25 
(middle  curve)  and  2.27cm  (upper  curve)  from  the  initial  interface.  The  initial  change 
in  temperature  is  zero,  but  the  middle  and  upper  curves  have  been  offset  by  5 and 
10  JiT,  respectively.  The  triangles  mark  the  beginning  and  the  end  of  the  current  pulse. 
An  thermoelectric  effects  are  included. 

Figure  3 : The  contributions  of  the  Peltier,  Thomson  and  Joule  effects  to  the  change  in 
temperature  as  a function  of  time  for  pure  Bi,  with  a current  pulse  of  80A/cm^  and 
duration  10s.  The  temperature  changes  are  measured  at  three  fixed  positions  in  the 
liquid,  0.05  (lower  curve),  1.25  (middle  curve)  and  2.27cm  (upper  curve)  from  the 
initial  interface.  The  initial  change  in  temperature  is  zero,  but  the  middle  and  upper 
curves  have  been  offset  by  5 and  10 A,  respectively.  The  triangles  mark  the  beginning 
and  the  end  of  the  current  pulse.  In  each  plot,  the  thermoelectric  effects  other  than 
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the  one  considered  are  set  to  zero. 


Figure  4 : The  contributions  of  the  Peltier,  Thomson  and  Joule  effects  to  the  interface 
velocity  as  a function  of  time  for  pure  Bi,  with  a current  pulse  of  80  A/cm ^ and  duration 
10s.  The  triangles  mark  the  beginning  and  the  end  of  the  current  pulse.  In  each  plot, 
the  thermoelectric  effects  other  than  the  one  considered  are  set  to  zero. 

Figure  5 : The  change  in  temperature  and  the  interface  velocity  as  a function  of  time 
for  pure  Bi,  with  a current  pulse  of  —SOA/cm^  and  duration  10s.  The  temperature 
changes  are  measured  at  three  fixed  positions  in  the  Hquid,  0.05  (lower  curve),  1.25 
(middle  curve)  and  2.27cm  (upper  curve)  from  the  initial  interface.  The  initial  change 
in  temperature  is  zero,  but  the  middle  and  upper  curves  have  been  offset  by  5 and 
10 K,  respectively.  The  triangles  mark  the  beginning  and  the  end  of  the  current  pulse. 
All  thermoelectric  effects  are  included. 

Figure  6 : The  contributions  of  the  Peltier,  Thomson  and  Joule  effects  to  the  change  in 
temperature  as  a function  of  time  for  pure  Bi,  with  a current  pulse  of  — 80A/cm^  and 
duration  10s.  The  temperature  changes  are  measured  at  three  fixed  positions  in  the 
liquid,  0.05  (lower  curve),  1.25  (middle  curve)  and  2.27cm  (upper  curve)  from  the 
initial  interface.  The  initial  change  in  temperature  is  zero,  but  the  middle  and  upper 
curves  have  been  offset  by  5 and  lOA^,  respectively.  The  triangles  mark  the  beginning 
and  the  end  of  the  current  pulse.  In  each  plot,  the  thermoelectric  effects  other  than 
the  one  considered  are  set  to  zero. 
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Figure  7 : The  contributions  of  the  Peltier,  Thomson  and  Joule  effects  to  the  interface 
velocity  as  a function  of  time  for  pure  Bi,  with  a current  pulse  of  —80 A/cm ^ and 
duration  10s.  The  triangles  mark  the  beginning  and  the  end  of  the  current  pulse.  In 
each  plot,  the  thermoelectric  effects  other  than  the  one  considered  are  set  to  zero. 

Figure  8 : The  change  in  temperature  and  the  interface  velocity  as  a function  of  time, 
and  the  solid  gallium  concentration  as  a function  of  position  for  an  alloy  of  germanium 
containing  0.1  at%  gallium.  A current  of  25  A/cm^  is  applied  at  t = 10s.  The 
temperature  changes  are  measured  at  three  fixed  positions  in  the  liquid,  0.02  (lower 
curve),  0.5  (middle  curve)  and  0.91cm  (upper  curve)  from  the  initial  interface.  The 
initial  change  in  temperature  is  zero,  but  the  middle  and  upper  curves  have  been  offset 
by  0.1  and  0.2Ar,  respectively. 

Figure  9 : The  change  in  temperature  and  the  interface  velocity  as  a function  of  time, 
and  the  sohd  gallium  concentration  as  a function  of  position  for  an  alloy  of  germanium 
containing  0.1  at%  galhum.  A current  of  — 25A/cm^  is  applied  at  t = 10s.  The 
temperature  changes  are  measured  at  three  fixed  positions  in  the  liquid,  0.02  (lower 
curve),  0.5  (middle  curve)  and  0.91cm  (upper  curve)  from  the  initial  interface.  The 
initial  change  in  temperature  is  zero,  but  the  middle  and  upper  curves  have  been  offset 
by  0.3  and  0.6K,  respectively. 

Figure  10:  The  change  in  temperature  as  a function  of  time  in  pure  germanium  and  in 
the  germanium  alloys  with  0.3  and  0.5  at%  gallium.  At  ^ = 10s,  a current  of  25  A/cm^ 
is  applied.  The  changes  in  temperature  are  calculated  at  distances  of  0.02  cm.  (lower 
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curve),  0.5  cm.  (middle  curve)  and  0.91  cm.  (upper  curve)  from  the  initial  interface 
position,  in  a moving  reference  system.  The  curves  are  offset  for  clarity. 

Figure  11:  The  change  in  temperature  as  a function  of  time  in  pure  germanium  and  in  the 
germanium  alloys  with  0.3  and  0.5  at%  galHum.  At  i = 10s,  a current  of  — 25 A/cm^ 
is  applied.  The  changes  in  temperature  are  calculated  at  distances  of  0.02  cm.  (lower 
curve),  0.5  cm.  (middle  curve)  and  0.91  cm.  (upper  curve)  from  the  initial  interface 
position,  in  a moving  reference  system.  The  curves  are  offset  for  clarity. 

Figure  12:  The  change  in  temperature  and  the  interface  velocity  as  a function  of  time, 
and  the  sohd  gallium  concentration  as  a function  of  position  for  an  alloy  of  germanium 
containing  0.1  at%  galhum.  A current  pulse  of  80A/cm^  is  applied  for  0.1  s.  The 
temperature  changes  are  measured  at  three  fixed  positions  in  the  hquid,  0.02  (lower 
curve),  0.5  (middle  curve)  and  0.91cm  (upper  curve)  from  the  initial  interface,  in  a 
moving  reference  system.  The  curves  are  offset  for  clarity. 

Figure  13:  The  change  in  temperature  and  the  interface  velocity  as  a function  of  time, 
and  the  sohd  gallium  concentration  as  a function  of  position  for  an  alloy  of  germanium 
containing  0.1  at%  galhum.  A current  pulse  of  — 80A/cm^  is  applied  for  0.1  s.  The 
temperature  changes  are  measured  at  three  fixed  positions  in  the  hquid,  0.02  (lower 
curve),  0.5  (middle  curve)  and  0.91cm  (upper  curve)  from  the  initial  interface,  in  a 
moving  reference  system.  The  curves  are  offset  for  clarity. 

Figure  14:  The  change  in  temperature  and  the  interface  velocity  as  a function  of  time,  and 
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the  solid  bismuth  concentration  as  a function  of  position  for  an  alloy  of  tin  containing 
0.78  at%  bismuth.  A current  pulse  of  50 A/cm^  is  applied  for  5.0  s.  The  tempera- 
ture changes  are  measured  at  three  fixed  positions  in  the  liquid,  0.03  (lower  curve). 
0.75  (middle  curve)  and  1.36cm  (upper  curve)  from  the  initial  interface,  in  a moving 
reference  system.  The  curves  are  offset  for  clarity. 

Figure  15:  The  change  in  temperature  and  the  interface  velocity  as  a function  of  time, 
and  the  sohd  bismuth  concentration  as  a function  of  position  for  an  alloy  of  tin  con- 
taining 0.78  at%  bismuth.  A current  pulse  of  — 50A/cm^  is  applied  for  5.0  s.  The 
temperature  changes  are  measured  at  three  fixed  positions  in  the  hquid,  0.03  (lower 
curve),  0.75  (middle  curve)  and  1.36cm  (upper  curve)  from  the  initial  interface,  in  a 
moving  reference  system.  The  curves  are  offset  for  clarity. 
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